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Abstract: We present global fits to the Large Volume Scenario (LVS) of string 
models using current indirect data. We use WMAP5 constraints on dark matter relic 
density, 6-pliysics and electroweak observables as well as direct search constraints. 
Such data can be adequately fit by LVS, with the best-fit point for /i > having 
= 13.6 for 8 degrees of freedom. The resulting constraints on parameter space are 
robust in that they do not depend much upon the prior, or upon whether one uses 
Bayesian or frequentist interpretations of the data. Sparticle masses are constrained 
to be well below the 1 TeV level, predicting early SUSY discovery at the LHC. 
We devise a method of quantifying which are the most important constraints. We 
find that the LEP2 Higgs mass constraint, the relic density of dark matter and the 
anomalous magnetic moment of the muon affect the fits to the strongest degree. 
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1. Introduction 

From the host of possibilities of Beyond the Standard Model physics the most well- 
studied and well-motivated is string theory. Stringy models naturally incorporate 
supersymmetry (SUSY), which can solve the problems of the instability of the Higgs 
mass (the hierarchy problem) and the nature of dark matter. Unfortunately the 
simplest extension of the Standard Model which includes SUSY, the so-called Mini- 
mal Supersymmetric Standard Model (MSSM) contains around 120 free parameters 
making a predictive analysis extremely difficult. While the effective dimensionality 
of the parameter space is reduced by constraints from flavour changing neutral cur- 
rents, it is still large enough to tax even the most powerful CPUs and techniques 
available today. It is to be hoped that some organising principle will be found at the 
LHC or a future collider to enable us to understand the relationships between these 
parameters, but in the meantime one popular way of dealing with this problem is 
by unifying the scalar mass terms to mo, the trilinear terms to A and the gaugino 
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masses to M at some high energy scale, usually taken to be the GUT scale ~ 10^^ 
GeV. With this pattern of soft SUSY breaking terms the MSSM is known as the 
Constrained MSSM (CMSSM) or minimal supergravity (mSUGRA). While univer- 
sality is a priori a very strong assumption, there are several string motivated models 
which predict such universality [1,2] and, more pragmatically, it renders the problem 
of performing a phenomenological analysis of the MSSM practical in computational 
terms. 

One such string theoretic model is the Large Volume Scenario (LVS) [2] in the 
limit of dilute fluxes. Originally discovered in the context of type IIB flux compact- 
iflcations, these models achieve 4D M = 1 broken supersymmetry with all moduli 
stabilised and exponentially large compactiflcation volume. All of these features are 
phenomenologically desirable. Having stabilised moduli means massless scalar par- 
ticles and non-realistic flfth forces are avoided. Exponentially large volume gives us 
confldence that the working in the supergravity limit is viable. Furthermore, the 
large volume V lowers the string scale and gravitino mass to 

Mp Mp 



From eq. ( |1 . 1|) it is clear that a volume V ~ 10^^ in string units will lead to TeV 
scale supersymmetry breaking as parametrised by the gravitino mass and to an in- 
termediate string scale rUg ~ 10^^ GeV. The beneflts of an intermediate string scale 
have previously been discussed in ref. [3] and include a natural solution to the strong 
CP problem, the correct scale of suppression of neutrino masses and gauge coupling 
uniflcation at m^, a possibility not usually considered in the context of mSUGRA. 
One further notable feature of the large volume models is that the flux superpo- 
tential Wq does not need to be flne tuned and is naturally of C(l), unlike the case 
of the well-known KKLT vacua [4] which require Wq ~ 10~^^. The Large Volume 
Scenario is therefore one of the most well studied and robust models which string 
theory presents. In this paper we shall confront the LVS with current indirect data 
from cosmology and particle physics, providing for the first time a global fit to the 
model where it is possible to trade goodness of fit between parameters, so that a bad 
fit to one observable may be compensated for by good fit to another. 

The connection between the high scale supergravity F-terms and the MSSM soft 
terms was established in [5] and demonstrated soft-term universality in the limit of 
dilute fiuxes. We have no knowledge of what the exact values of the high scale bound- 
ary conditions are, and a systematic exploration of the parameter space has yet to be 
performed. Such an exploration would tell us what the viable regions of parameter 
space are and equally as important, what regions have already been ruled out. We 
also wish to take into account uncertainties and errors in our knowledge and pre- 
dictions about standard model parameters. The usual procedure of fixed grid scans 
through parameter space is therefore not useful, as well as being computationally in- 
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tensive since the number of points scanned and hence the time taken is proportional 
to , where is the number of important free parameters the model has and k is 
the desired number of points along each dimension. In standard mSUGRA this is 
taken to be 8 and in our model it is 6. A more efficient method that allows us to 
take errors and uncertainties into account is the use of Monte Carlo Markov chains 
(MCMC) in a Bayesian statistical formahsm, as was first considered in ref. [7] and 
further developed in refs. [7-10]. 

What we are interested in is the posterior probability distribution function (pdf ) 
p(m|data), the probability of a point m in parameter space being "correct" given 
some Standard Model data such as masses and other observables. Unfortunately is 
difficult to calculate: given a value of some branching ratio it is arduous to invert this 
to find which model parameters are consistent with it. However, it is quite straight- 
forward to obtain the likelihood distribution p(data|m) which is the probability of 
obtaining some particular SM observables given a point in parameter space. Once 
we have the spectrum of the model at that point it is trivial to calculate things like 
branching ratios. The connection between the likelihood distribution and the poste- 
rior pdf has been provided for us by the Reverend Thomas Bayes and his eponymous 
theorem [11] 

p(mi|data) p(data|mi)p(mi) 
p(m2|data) p(data|m2)p(m2) 

where p{mi) is known as the prior (pdf) or simply prior. It encodes our previous 
beliefs or uncertainties about a particular point in parameter space. In the case of 
flat priors where p(mi) = p{m2) the likelihood distribution is proportional to the 
posterior pdf. We wish then to numerically construct the likelihood distribution. 
To do this we use Monte Carlo Markov chains. A Markov chain is best described 
as a string of points which sample from some continuous distribution. These have 
the benefit that the run time depends only linearly on the number of dimensions in 
parameter space, in contrast with the power law behaviour for grid scans. For more 
details on the implementation of the MCMC method we use in this paper see ref. [7] . 
To construct the hkelihood we will use a variety of data from cosmology, electroweak 
precision observables, ft-physics and current sparticle mass limits, including for the 
first time in a Bayesian context BR(B rv), the meson mass splitting AMb.. and 
the isospin asymmetry Aq-. We also make use of the new WMAP5 dataset and the 
most recent measured value of the top quark mass rrit from the Tevatron. 

Previous phenomenological studies of the LVS are refs. [5,6] where random sam- 
ple spectra were generated which were then used to place bounds on the viable regions 
of parameter space and investigate LHC coUider observables and signatures of the 
model. This paper does not address the issue of collider observables, but extends 
the above cited work to accurately fit and sample the parameter space of the model 
using currently available indirect constraints. 
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Recently rcf. [12] have performed a analysis of the CMSSM and minimal 
GMSB and AMSB scenarios using a subset of electroweak and 6-physics observables. 
Their fits consist of random scans of parameter space using 10^ points. Aside from 
the fact by incorporating variations in the SM parameters and, importantly, the dark 
matter rehc density which is known to be the strongest constraint on the CMSSM 
parameter space, our MCMC approach also allows for better sampling and hence 
statistically more stable results. On the other hand, we do not do multi-model 
hypothesis testing, but instead perform a hypothesis test on the sign of fj,, the bi- 
linear parameter of the Higgs potential whose magnitude is accurately known by 
fitting to Mz- 

In the next section we discuss the origin of the soft SUSY breaking terms in the 
Large Volume Scenario and some caveats regarding some approximations we make. 
In Section 3 we present our suite of observables with which we will constrain the 
model and describe how we construct the likelihood, the issue of priors and how we 
test for convergence of our Markov chains. Section 4 presents our fits to the hkelihood 
distribution and posterior pdfs for some sparticles and the dark matter relic density. 
We also discuss channels of relic density annihilation, best-fit points for the model 
and present a new variable constructed to make quantitative statements about which 
observables are constraining the form of the likelihood distribution the most. Finally, 
we entertain the possibility that /x < 0. Section 5 describes a frequentist's approach 
using profile likelihoods which we extract from our Markov chains and identifies 
some interesting 'volume effects'. We conclude by recapitulating our main points 
and presenting some possible directions for future research. 

2. Soft terms in the LVS 

The Large Volume Scenario [2] is one of the most successful paradigms in string 
theory which achieves realistic low energy physics with stabilisation of all moduli 
fields. From within this scenario one can obtain TeV scale supersymmetry (SUSY) 
breaking, inflation [13], QCD axions [14], and the correct scale of neutrino masses [15]. 
For a comprehensive review see [16]. The scenario arises by considering generic 
quantum corrections to the string action in KKLT style compactifications [4] on a 
Calabi-Yau manifold with a "Swiss-cheese" style geometry, and is immune to the 
usual flne-tuning problems of KKLT compactiflcations. There is a danger that once 
we consider some quantum corrections it is inconsistent to ignore all the others: this 
is the Dine-Seiberg problem [17]. However, initial studies suggest that the models 
are robust against further corrections [18,19]. 

When compactified on a Calabi-Yau orientifold type IIB string theory can be 
described at low energies by an effective J\f — 1 4D supergravity theory. In this paper 
we assume that this theory has the matter content of the Minimal Supersymmetric 
Standard Model (MSSM). The Lagrangian is then determined by a Kahler potential 
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/C, a superpotential W and gauge kinetic functions fa, which we may expand in terms 
of moduh fields $, matter fields C" and the two MSSM higgs fields Hi and H2 as 
follows: 

W = W{^) + ii{^)HiH2 + ^ra;3^($)C"C"^CT + . . . , (2.1) 
K = + $)C"C^ + (Z($, $)//ii72 + /i.e.) + . . . , (2.2) 

/a = /a($). (2.3) 

We will review the origin of the MSSM soft terms but will gloss over the more 
technical aspects of the derivation which may be found in [5]. We turn first to the 
gauge kinetic functions, gauge couplings and gaugino masses. 

2.1 Gauge Couplings 

To be considered truly complete, any stringy model should have a sector which 
contains the Standard Model. We assume that such a sector can be found where the 
Standard Model will be supported on magnetised D7 branes. In order that the SM 
gauge groups not be too weakly coTiplcd, these branes must wrap a small 4-cycle Ts 
in the Calabi-Yau. The gauge kinetic functions may be computed by dimensional 
reduction of the DBI action whereupon one obtains 

/a = ^ + K{F)S, (2.4) 

where S is the axio-dilaton field, is the Kahler modulus of the small cycle and ha 
is a topological function of the fluxes present on the D7 brane. These functions are 
in general dependent on the explicit brane configuration used to realise the Standard 
Model are and currently unknown for realistic scenarios outside of toroidal orien- 
tifolds [20] If the cycle size is increased, the fluxes become diluted and the gauge 
couplings become independent of the fluxes. In this dilute flux limit we may then 
write 

/Sf7(3) = 

fsui2) = i, (2.5) 

Ml) = ^Yii, 

where ky is a generally model dependent normalisation for the U{1) gauge field, 
which we regard as unknown. In the dilute fiux limit this parameter will not affect 
the physics. 

The gaugino masses are 

^" - -2-1^' P-^' 

where are the supergravity moduli F-terms 

pm ^ f,K/2j;.mnj^_^ (2.7) 
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which quantify the supersymmetry breaking. In the dilute flux limit this gives 



Ml = M2 = M3 = — = M, 
2t. 



(2.8) 



where Mj = Msu{i), so that the gaugino masses are universal at the compactification 
scale, which in these models is at the intermediate scale ~ lO^^GeV. 

2.2 Soft Terms 

To derive the rest of the soft terms we use following standard expressions for the 
scalar masses, trilinear A-terms and the S-term [21]: 



< = + Vo) - F'^F'^dndmhgK^, 

Aaf^-y = F"^ K^dm log Y^p^ - dm \og{k^kpk^ 



Bfi = {KH.KH.r^'H e^/^ fF"^ \km + dm log/i - d\og{kH,kH, 



(2.9) 
(2.10) 



"^3/2 



+ (2m^ ,2 + V^)Z - m^/^F'^dmZ + ms/sF'" d^Z - Zd^ logikn^ku. 



dmdnZ - {drnZ)dnlog{KH,KH2] 



(2.11) 



where K^p = KaS^p ( without summing over a) is the matter field metric. Using 



the metric for chiral matter fields derived in ref. [22] 



" V2/3' 

leads to the soft terms in the dilute flux limit [5] 

M„ = M 
m„ = a/am 

AalB'y = -3AM 

B = -(A + 1)M 



(2.12) 



(2.13) 



where A is the modular weight of the matter fields with respect to the small cycles. 
In the minimal case where all branes are wrapping the same cycle it was shown in 



ref. [22] that A = 1/3 so that eq. p.l3| reduces to 



M, = 


M 


m„ = 


M 






-M 


n 


AM 
3 



(2.14) 



which interestingly reproduces the form of soft terms in the dilaton-dominated sce- 
nario of heterotic string models. It is this minimal case which we will consider in 
this paper. 
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It was demonstrated in [5] that it is impossible to implement the 5-term condi- 
tion for /i > 0. The condition can be satisfied for /i < but only in a region of low 
tan (3 which leads to a Higgs mass below the lower bound from LEP. We must then 
assume that there exists a means to generate the correct 5-term, for example an 
NMSSM-style coupling aNHiH2, where is a gauge invariant scalar which obtains 
a vev. There also might exist vector-like matter between the string scale and the 
TeV scale which could alter the RG equations and the low-scale soft terms. Perhaps 
the most important caveat is that we are using universal gaugino masses (|2.8|) at 
the intermediate scale. It is well known that the standard model gauge couplings 
with MSSM field content unify at the GUT scale. It is clear from eq. ( p.4|) that 
the required non- universality is provided in the LVS by the fiuxes. We can estimate 
the magnitude of the effect of the fiuxes by running the Standard Model SU (2) and 
SU (3) couplings to the intermediate scale and noting the non-universality, obtaining 
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Ms 
M2 



1.37. (2.15) 



We adopt a compromise position of leaving the gaugino masses as universal at the 
string scale, but allowing the gauge couplings to differ by the amount above. In hght 
of this discussion all results in this paper should be understood as being to leading 
order in the dilute flux approximation (p.5|). 



3. Observables and the Likelihood 



In calculating the likelihood we 
follow refs. [7] and [23]. We vary the 
six inputs shown in Table |l[ Several 
runs with varying starting points in 
parameter space were performed with 
an expanded parameter range 60 < 
mo < 2000GeV and 2 < tan/5 < 62 
which showed that consistent solu- 
tions of the RG equations only ex- 



mSUGRA parameter 


range 


mo 
tan (3 


60 GeV to 1.5 TeV 
2 to 30 


SM parameter 


constraint 


mhirrih)^^ 


127.918±0.018 
0.1172±0.002 
4.20±0.07 GeV 
172.6±1.4 GeV 



isted within the more limited region Table 1: Input parameters 

shown. To increase the efficiency of 

our simulations we restrict these parameters accordingly. The four SM inputs are 
constrained to lie within 4a of their central values, l/a'^^ and m;,(mf,) are taken 
from ref. [24], a^^ from ref. [25] and rrit from the most recent combined Tevatron 
analysis [26]. Given how small the experimental errors are we flx the muon decay 
constant to be 1.16637 x 10^^ GeV~^ and the mass of the Z vector boson to be 
91.1876GeV [24]. To calculate the spectrum of the MSSM we use a modifled version 
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of SOFTSUSY 2.0.17 [27]. The W mass and sin^ are obtained with a code based 
on refs. [39,40]. 
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67.7 


m-g 


195 




76 




88 
43.1 




86.4 




91 




250 



Table 2: Lower bounds applied to sparticle mass predictions (in GeV). 

If a generated spectrum contains a sparticle whose mass hes below the 95% 
lower bounds listed in Table 2 [9] it is assigned zero likelihood. Similarly, if a point 
contains tachyonic sparticles or does not break electroweak symmetry in the correct 
way it is assigned zero likelihood. If a point has survived this far it is passed via 
the SUSY Les Houches Accord (SLHA) [28] to micrOMEGAS 1.3.6 [29] to calculate 
the dark matter relic density, the branching ratio of the rare decay BR{Bs —>■ 
and the anomalous magnetic moment of the muon, {g — 2)^, and to the most up-to- 
date version of Superlso [30] to calculate the the isospin asymmetry of the decay 
B K*'y and the branching ratio BR{b 57). 

3.1 Observables 

To construct the likelihood we use a collection of observables from across electroweak 
physics, cosmology and B-physics. 

We use the newly released 5- year WMAP 
data [31] to constrain the cold dark matter 
relic density. 



50 



40 



QnMh^ = 0.1143 ±0.0034 



(3.1) 



As has been recently pointed out in ref. [32] 
small and currently undetectable changes in 
the expansion of the universe before Big- 
Bang Nucleosynthesis (BBN) can lead to in- 
creases (but not decreases) in the relic den- 
sity by factors up to 10^. We must also con- 
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our constraint 


pure WMAPE 




/ \ 





0.05 0.1 



0.15 0.2 



Figure 1: Constraints for the dark mat- 
ter relic density 

sider the possibility that the neutralino does not constitute the entire relic density: 
there may be other components such as axions or dark fluids. More recently, ref. [33] 
have discussed the effects of massive right-handed neutrinos on the mSUGRA relic 
density. Contrary to common lore, the neutrino Yukawa couplings can increase the 
relic density by more than an order of magnitude. We therefore adopt the suggestion 
of [32] and given a prediction of the relic density u we take the likelihood to be 



1 



c + cr 



{nh^ < 0.1143), 



-(c-tj)2/2(T2 



c + a 



{Vth^ > 0.1143). 



(3.2) 



- 8 - 



In this way points below the WMAP bound are assigned an equal likelihood, whereas 
those above are subject to a Gaussian likelihood centred on c = 0.1143 with standard 
deviation a = 0.02, which represents theoretical error in the prediction of the relic 
density. We show this in Fig. |l|, along with the constraint we would obtain from pure 
WMAP5. 

The anomalous magnetic moment of the muon gains contributions from QED, 
hadronic vacuum polarisations and light-by-light processes. Evaluating the vacuum 
polarisation with e^e~ and r data gives different results. Since up to date e~^e~ 
results are in line with earlier results, and those for the r data are not, we restrict 
ourselves to the e+e~ data and obtain [34] 

a^^^ = (11659178.5 ± 6.1) x 10"^°. (3.3) 

When compared with the experimental result 

a^^P = (11659208.0 ± 6.3) x 10"^° (3.4) 

there exists a discrepancy 

^ (^_^2V ^ ^ ^exp _ ^SM ^ (29.5 ± 8.8) X 10-^° (3.5) 

at the 3.4cr level, which may be an indication of new physics. We evaluate this at 
one-loop with micrOMEGAS, and then add in the logarithmic piece of the QED 2- 
loop calculation and the 2- loop stop-higgs and chargino-stop/bottom contributions 
[36,37]. 

To the experimentally measured mass of the W vector boson and the effective 
weak leptonic mixing angle [38] we add the SM and MSSM theory errors detailed 
in [39,40] to obtain 

Mw = 80.398 ± 0.027 GeV, sin^ 9^ = 0.23149 ± 0.00017. (3.6) 

We utilise the full MSSM one-loop contribution, SUSY corrections of order 0{aas) 
and 0{a^fj), as well as all relevant SM-like terms beyond one- loop order (see refs. 
[39,40] for details). 

We also use the LEP2 Standard Model Higgs mass bound that rrih < 114.4 GeV 
at 95% confidence level. We smear this by 3 GeV to represent the theoretical uncer- 
tainty in the SOFTSUSY prediction for rrih, as described in ref. [23]. 

Flavour-changing neutral currents are one of the areas most sensitive to new 
physics, particularly in the large tan/3 regime. With the accumulation of precision 
results from BABAR and Belle over recent years these observables have gained con- 
siderable discriminatory power. 

The rare branching ratio BR{b s'j) has been measured as [41] BR{b — > sj) = 
(3.55 ± 0.26) X 10"^ We add in the SM and MSSM uncertainties as in ref. [42] to 
obtain 

BR{b S7) = (3.55 ± 0.72) x 10"^ (3.7) 
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For the branching ratio BR{Bs jJ,^ jJ. ) the most recent upper bound from the 
Tevatron is [43] 

BR{Bs i^^iJ.-) < 5.8 X 10"^ (3.8) 

at 95% confidence level. 

The current HFAG average of the branching ratio for the process B^ tv is [44] 

BR^'^PiBu Tv) = (1.41 ± 0.43) x 10"! (3.9) 

The Standard Model prediction of this branching ratio depends on whether one deter- 
mines the CKM matrix element \ Vui,\ inclusive or exclusive of semileptonic decays [45]. 
We statistically average over these two values and get BR^^ = (1.12 ± 0.25) x 10"*^, 
so that 

1.259 ±0.378. (3.10) 

The MSSM contribution to this branching ratio is dominated by charged-Higgs con- 
tributions and to leading order in tan {3 gives 



R 



Btv 



1 - 



m% \ tan^/3 



?77.J^± / (1 + Co tan (5) 



(3.11) 



where eo is an effective coupling which takes into account the non-holomorphic cor- 
rection to the down-type Yukawa coupling induced by gluino exchange [46] and is 
given by 

2as^, (Ml M 

to — — 

where 



H2{x, y) = r + r, (3.13) 

{l-x){x-y) {l-y){y-x) 

The mass splitting of the Bs meson has been measured by CDF to be [48] 

A'^^Ms^ = 17.77 ± 0.12ps-\ (3.14) 

while the UTFit evaluation of the standard model estimate is A^^m, = 20.9 ± 
2.6ps~^ [49]. This gives us 

^AM.. = = 0.85 ± 0.12. (3.15) 

The dominant MSSM contribution to AMg^ comes from neutral Higgs particles in 
double-penguin diagrams [47] and is given by 

Rambs ^ '^-n^b{mb)ms{mb) 647rsm 6'gff (eyA^ tan /5) 



aemMlSoimj/ml^) [i + (e^ + eY\j)tan(3f [1 + eotan/5]^' 

(3.16) 
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where 

Af (Ml Ml , 
6y = -T^H^ M^, ^ (3.17) 



and Co is the same as in eq. ( |3.12| ) and 5*0 is a Wilson coefficient which can be found 



in ref. [50]. This does not take into account the charged Higgs and chargino box 
diagrams, however these provide only small contributions to the total splitting over 
the majority of parameter space. 

The Standard Model predictions for the two above observables both depend in 
some way on the mass of the bottom quark. Since this is also something that we will 
be predicting and fitting to in our model, we must consider the possible correlations 
between the Standard Model and SUSY contributions. The main parameters which 
could depend on are the B-meson mass and the mesonic decay constant /s. 
We take the from experiment to be 5.279 GeV [24], and so is not correlated with 
any particular value of m^. The decay constant is calculated on the lattice and 
the associated value of mf, is not an input but is derived by fixing a scale after which 
predictions may be made. The value of m;, obtained in this way is consistent with 
all values that we will consider in this paper, see ref. [51]. 

Our final observable is the isospin asymmetry from the exclusive process B 
K*7, defined as 

. _ r(E°^ir*°7)-r(i?-^ir*-7) 

Data from BABAR [52] and Belle [53] have constrained this to be 

-0.018 < Ao_ < 0.093 

at 95% confidence level. We convert this to a Gaussian with central value and 
standard deviation 

Ao_ = 0.0375 ± 0.0289. (3.19) 

3.2 The Likelihood 

A prediction of one of the above observables (except QoAih'^, BR{Bs Ai^/^ ) 
and rrih which we treat separately) is assigned the log likelihood 

jc j-pi)'^ _ 1 
2s^ 2 

where Cj is the measured central value and Sj the standard deviation. The relic 
density likelihood is treated as in eq. (|3.1|) . The likelihood for the branching ratio 
of Bs —>■ fJ'^fi" is calculated using the predicted value from micrOMEGAS and CDF 
Tevatron Run II data [54]. For the Higgs we use a parametrisation of the search 
likelihood from LEP2, as used in ref. [55]. We then calculate the combined likelihood 
for all observables 

ln/:*°* = ^InA, (3.21) 



In A = - -ln(27r) -ln(s,) (3.20) 
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equivalent to the assumption that our observables form a set of independent quanti- 
ties. 



3.3 Priors 

In a Bayesian framework our ignorance is quantified through the posterior probability 
density function. Previously [7, 9] , this has been set to be 

p(mo, tan/?, .[data) = p(data|mo, tan/3, - ^) , (3.22) 

p(^Qata J 

where s are some Standard Model inputs and j9(data|motan/5, s) is the likelihood. If 
we desire the posterior pdf for a particular parameter we marginalise over (i.e. inte- 
grate out) all other parameters. The natural measure which we use for marginalising 
comes from (|3.22 ) so that if we wish to know the posterior pdf for mo, for example, 
we calculate 

p(mo|data) = J rftan/^rfs p{mQ, tan P, s\data). (3.23) 

It was observed in [8] that this does not reflect the fact that tan/3 is a parameter 
derived from the more fundamental parameters B and /i. What we really desire is 

p(mo|data) = J djjdBds p(mo, tan/3, s|data)5(M^ — M^^^) 

dtanPds p(mo, tan/5, s|data)r(i?, yU, tan/3) , (3.24) 



where we have fixed Mz at its experimental value since the experimental error is so 
small and r{B, /i, tan f3) is a Jacobian factor. To determine r we follow ref. [8] and use 
the relations between tan/5, Mz, B and /i derived from the electroweak symmetry 
breaking conditions [56] 

sin 2/? , 2 , - 2 , o 2\ 
fiB = [niH^ + m^^ + 2/i ) , 



m 



- mjj^ tan^ /? M| 



^2 = .7 ' - (3.25) 



tan^ /5 - 1 2 



to obtain 



r(5,/i,tan/5) = Mz 



B tan^ /5 - 1 



(3.26) 



/itan /? tan^ /5 + 1 

which will from now be referred to as the radiative electroweak symmetry breaking or 
REWSB prior. Viewed as a measure on parameter space, the REWSB prior assigns 
a higher weight to points with lower tan /5 and [i. 

3.4 Convergence 
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To accurately sample the posterior 
probability density we ran 10 indepen- 
dent MCMCs of length 10^ each. We dis- 
carded the initial 4000 steps as "burn-in" 
for the MCMCs. We use the Gelman- 
Rubin R statistic [7,57] to check for con- 
vergence. In this test convergence is in- 
dicated by the value r < 1.05. With 
this measure all of our runs converge in 
less than 10,000 steps, and reach final 
R- values of 1.002 after 10^ steps. We 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 O.S 

No. of points /lO^ 



Figure 2: No. of points plotted against R, 

the Gelman-Rubin convergence statistic, 
illustrate this in Fig. ^ by plotting the 

Gelman-Rubin statistic against the number of steps taken. Convergence is achieved 
here very rapidly, in around 10^ steps. This should be contrasted with the mSUGRA 
case in ref. [7], where it took 6 x 10^ steps to converge. This is because the mSUGRA 
likelihood distribution is a complicated multimodal distribution spread out over a 
large region of parameter space, unlike the relatively compact LVS likelihood distri- 
bution we will see below. The overall efficiency of our simulation is 44% (43%) for 
/i > and 23.9(22.8)% for < with flat (REWSB) priors respectively. Throughout 
this paper when binning we use 75 x 75 bins, and all 2D plots are normalised to the 
maximum likelihood bin. In all following ID plots (except in Section 5) the vertical 
axis is the posterior probability per bin. 



4. Likelihood Fits 

Here we present the results of our analysis for values of /i > 0, with flat and REWSB 
priors. In Fig. ^(a) and (b) we present the posterior pdf marginalised to the mo-tan/? 
plane. We also show 68 and 95% confidence limit contour lines. Since M1/2 and A 



are determined linearly in terms of mo according to eq. ( |2.14]) , we do not present 
these pdfs. The sharp cutoff at around tan (3 ~ 17 is due to points to the right of 
the cutoff having the stau as the LSP, which have been rejected with zero likelihood. 
Fig. ^(c) shows the ID likelihood distribution of the variable mo. The majority of 
the likelihood is located between 200 and 400GeV for both sets of priors, with the 
upper bounds at 95% c.l. being 389.2GeV(616.2GeV) for flat (REWSB) priors. The 
REWSB priors have a fatter tail to higher mg due to the boomerang shape of the 
posterior pdf, and the REWSB priors favouring lower tan/3. Fig. |^(d) shows the 
likelihood distribution for tan (3. Solutions to all constraints exist only in the range 
2.5 < tan/3 < 20, and the 95% upper bounds are 17.1(16.0) for flat (REWSB) priors. 
Although tan (3 is difficult to determine experimentally, if measurements showed that 
we live in a high tan/3 region (as favoured by mSUGRA for example [7]) this would 
discriminate against the Large Volume Scenario. 
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Figure 3: Marginalised posterior pdfs in the m-o-tan /? plane for /i > and (a) flat pri- 
ors and (b) REWSB priors with 68% and 95% c.l. contours shown, and the binned ID 
distributions for (c) uiq and (d) tan/3 also for ^ > 0. 



In Table ^ we show the details of the best fit points for the flat and REWSB 
priors. For mo and tan (3 we show the values of these parameters, and for all other 
quantities we give the that quantity makes to the total of the fit. The third 
column presents our estimate of the uncertainty we have in our estimate of the 
absolute minimum of the minimum of all examples for the case of flat priors. To 
compute this we treat the best-fit point of each of the ten chains as an independent 
estimation of the minimum for the observables. From these ten points we calculate 
the standard deviation of the minimum value, which we take as the uncertainty 
in our estimates. Given an infinite run time we would expect the flat and REWSB 
priors to converge to the same best-fit point. The fact that these are different is 
merely indicative of the finite length of our Markov chains. The neutralino makes up 
around two-thirds of the dark matter relic density and the lightest CP-even Higgs 
mass is in both cases around 115 GeV, just above the lower bound set by LEP2. 
Both sets of priors slightly underpredict by about la. AMb^ varies by less 
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-0.48 
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Mw 


0.94 


1.12 


0.04 


{g - 2)m 


5.48 


9.54 


0.20 


Ao_ 


3.55 


2.31 


0.16 




1.56 


1.56 





BR{Bu Tu) 


0.57 


0.48 


0.01 






(total) 


13.6 


15.3 


0.83 



Table 3: Best fit points for /i > and statistical pull of observables. mg and tan/3 are the 
parameter space values we generate. All other numbers are x' values. The third column 
gives the statistical error on our estimates of the x' values for the flat prior best-fit point, 
as described in the text. 



than one percent over the entire parameter space, and Rambs never less than 0.99. 
Similarly, the branching ratio BR{B tu) does not exhibit much variation, and 
Rbtu is always above 0.87, inside the la bounds we derived in Section 2. This is 
in agreement with what was found in ref. [58]: that these two observables do not 
impose large changes to the fitted parameters. The fiat priors have an overall slightly 
smaller x' value, indicating a better fit to all observables, and so for the remainder 
of this paper we take this to be the "true" best-fit point. The REWSB priors have 
an overall better fit to the B-observables: a x' value of 4.6 versus 7.1 for the flat 
priors. However, due to the value of tan/3 in the REWSB case being just under half 
that of the flat case, the supersymmetric contributions to the anomalous magnetic 
moment of the muon are very small for REWSB, which leads to the best-fit point 
of the fiat priors having a smaller overall x' value. The reason that the relic density 
has no x' associated with it is because the value of the best-fit point falls below the 
WMAP central value, and does therefore not incur a likelihood penalty as discussed 
in Section 3.1. 

We may perform a simple hypothesis test by computing the P-value associated 
with the x' of the best-fit point. To do this requires that we know the number 
of statistical degrees of freedom of our fits. We argue that this is eight: we have 
fourteen observables defined above (if one includes m^, m^, and a~^) and six 
model parameters: mo, tan/5 and the four SM observables just mentioned which are 
varied in our MCMC runs. The P-value of the best-fit point is then Pfiat = 0.093. 
Taking a standard significance level of 0.05, we may then say that the LVS is indeed 
consistent with the available data. 

4.1 Dark matter 

In Fig. ^ we show the ID likelihood distribution for the dark matter relic density. 
The most striking feature for both sets of priors is the large spike in likelihood near 
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the origin. To elucidate what is responsible for this we turn to discussing channels 
of relic density depletion. 

We would like to assign probabilities 
to the possible relic density depletion pro- 
cesses. To this end we follow ref. [7]: the 
stau co-annihilation region is where m^o lies 
within 10% of m^^, the h^/A^/Z pole re- 
gion is where 2m^o is within ten percent of 
m/jo/m^o/Z, respectively, and the stop co- 
annihilation is situated where m'j' is within 
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since this channel is particu- 
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Figure 4: 1-D posterior distribution for 
the relic density for flat and REWSB pri- 
ors for > 0. 



larly efficient. Since the lightest neutralino 
is always more massive than 120 GeV the 
/i°-pole and Z-pole regions are completely 
inaccessible throughout the parameter space 
Similarly, there is a negligible amount of 
stop CO- annihilation. However for both sets 
of priors the probability is higher than 99.5% 

that we are in both stau co-annihilation and A-pole regions. This is in sharp contrast 
with mSUGRA, where the pole and stau co-annihilation regions do not overlap 
and where the A^ pole region is only found at high tan (3. 

To further investigate the spike we filter our chain of points, keeping only those 
with QDjiih"^ < 5 X 10~^. To obtain such a small relic density requires (at least) one of 
the annihilation channels to become extremely efficient, which should occur when the 
masses of some of the neutralino and some of the particles above become degenerate. 



We therefore plot in Fig. H(a) the mass differences m 
filtered spike region, and for both sets of priors. 

The stau-neutralino mass splitting is 
always less than 25 GeV, peaking at just 
less than 20 GeV. This peak is more pro- 
nounced with the REWSB prior and is 
associated with the spike region. How- 
ever, the spike pdf does not extend down 
to the degenerate region of the plot and 
stau co-annihilation does not significantly 
contribute to the relic density depletion in 
the spike. 
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Table 4: Top five relic depletion channels 
for the best-fit points in Table |3| for /x > 0. 



The 



In Fig. ^(b) we show the distribution of the mass difference tua — 2m^o. 
maximum of the plot occurs near —20 GeV for the flat and REWSB priors, and for 
the flltered spike region the peak is nearly exactly at zero. The implication of this 
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Figure 5: Pdfs for the mass splitting between (a) Xi ^ (b) between tua and 
2m^o with flat and REWSB priors, and for the filtered spike region discussed in the text 
for fi > 0. 

is that the spike is associated with a region where ttia ~ 2m^o. There still remains 
the question as to whether the spike is a region of high likelihood, or a large region 
of average or low likelihood. As we shall show in section |^, it is the latter that is the 
answer, and that the spike is an example of a so-called volume effect. 

Table ^ shows the top 5 relic depletion 
channels for the best fit points above, as cal- 
culated by micrOMEGAS. The channels XX ~^ 
bb, ti and rf are all s-channel interactions 
mediated by the Higgs. These channels 
are responsible for 68%(82%) of the relic den- 
sity depletion for the fiat (REWSB) best-fit 
points. This is in agreement with what we 
would expect from Fig. |[ Subdominant pro- 
cesses include stau co-annihilation, selectron 
co-annihilation and ff tt. 

Although we have purposefully left open 

the door open for exotic dark matter in our -r,. a t3 ^ ■ Af • v ^ + 
^ b igure d: Foster lor pdt marginalised to 

choice of rehc density constraint, the scenario ^he mo-tan /? plane with WMAP5 dark 
that is both minimal and most discussed in j^^tter constraint and /i > 0. 
the literature is where the relic density is en- 
tirely neutralino. To allow comparison with this case we reweight our Markov chains 
so that likelihood function is a symmetric Gaussian centred on 0.1143 with stan- 
dard deviation 0.02, which is the pure WMAP5 constraint shown in Fig. |I|, and then 
re-analyse all our data. Fig. ^ shows the resulting likelihood distribution in the rrio- 
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tan P plane. As expected, the distribution is much more constrained when compared 
with Fig. ^, with most of the loss coming from the low mo region. Points in the 
low rriQ region have increased dark matter annihilation cross-sections due to lower 
sparticle masses and so the relic density tends to be low, and it was this area which 
was responsible for the spike in Fig. H. 
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Figure 7: ID Mass distributions for (a) the gluino (b) squark g^, (c) selection (d) 
the CP odd higgs , (e) the lightest neutralino Xi, (f) the stop t and (g) the lightest higgs 
h for flat and REWSB priors with fi > 0, and with profile histograms. Profile likelihoods 
are discussed in Sect. ^, and have been rescaled to aid visual comparison. 

Fig. 1^ shows the likelihood distributions for the masses of (a) the gluino, (b) the 
left-handed squark, (c) the right-handed selectron, (d) the CP odd Higgs A^, (e) the 
lightest neutralino, (f) the stop and (g) the lightest CP even higgs h for fi > with 
flat and REWSB priors, and for the profile likelihoods discussed in Sect. ^. Table 
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1^ shows the 95% confidence hmit upper bounds on these masses. It is notable that 
the 95% upper bounds are almost all less than 1.5 TeV and therefore lie well within 
the reach of the LHC within the next few years. The higgs mass rrih is constrained 
to lie below 120 GeV, implying a late discovery of this particle at LHC. The priors 
have only a small effect on the shape of the distribution, and our predictions can 
therefore be considered quite robust. 
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Table 5: 95% c.l. upper 
bounds for sparticle masses for 
/i > 0. All figures are in GeV. 



the tail out to higher masses. This is because the 
REWSB prior favours lower tan j3 which pushes the 
favoured region of parameter space up the tail of the 
"boomerang" in Fig. ^(b) to higher (and there- 
fore higher mi/2), thereby increasing the likelihood 
of higher sparticle masses. Also of note is that the 
peaks of all three distributions occur very close to 
one another. This independence of the posterior pdfs 
from the priors and agreement with the profile like- 
lihood indicates that there are enough observables 
constraining our model to overcome whatever prior 

beliefs we might have. The predicted masses are therefore indicative of the LVS and 
will not change with addition of more data. 

It is also interesting to quantify which observables are constraining the likelihood 
the most. An observable which is essentially constant over the parameter space 
will not constrain the posterior pdf very much, while one which exhibits significant 
variation over a range of a few standard deviations from the experimental central 
value would contribute a large to the likelihood in some regions but very little 
in others. To see how an individual observable is constraining the posterior pdf for 
the model variables mo and tan (3 we consider two cases: when it is used in the 
construction of the likelihood, and when it is omitted. To correctly take into account 
correlations between observables we consider the posterior pdf of mo jointly with 
tan /3, whose volume we normalise to one. We then calculate the integrated posterior 
difference of tan /3 jointly with mo which we call the "moulding power" 



Mp = 1/2 J d tan/5(imo |p(all data|mo, tan/3) — p(all but one data|mo,tan/5)| . 

(4.1) 

It is this quantity which we use as a measure of the effect an observable has 
on the likelihood, and the fact that each of the ten chains provides a statistically 
independent determination of Mp. Table ^ shows the estimates for Mp obtained by 
this procedure. As expected BR{Bs — > fifi) and AMp^ do not constrain the form 
of the likelihood. Similarly, Mw and sin^ are effectively constant on parameter 
space. The main constraining observables are the Higgs mass rrih and the relic 



- 19 - 



density QoAih^, with the anomalous magnetic moment {g — 2)^ being the next most 
constraining. 

We know that in mSUGRA that 
the unconstrained rehc density can reach 
values as high as 100 [7], and with the 
Gaussian constraint we use, accurate sam- 
pling of the high relic density region is 
not possible for distributions with nar- 
row allowed regions and long tails. We 
found that the tail will was not accu- 
rately sampled, leading to bad statistics 
after reweighting resulting in a large value 
standard deviation of Mp. We there- 
fore ran 10 chains of length 50000, omit- 
ting the relic density from the construc- 
tion of the likelihood and used these Table 6: Moulding power of individual ob- 
chains to calculate p(all data but nDMh^\mSjM''p)^ described in the text, for ^l > 0. 
In sharp contrast with mSUGRA we found 

no points where the relic density was higher than 1.5. For comparative purposes we 
have also calculated the constraint value for the relic density with the pure WMAP5 
constraint in Table ^. 

Finally, we comment on the "golden channel" decay chain Ql ^ X2 ^ ^ 
Xi- This important chain can give constraints on the mass spectrum [59] and even 
information about sparticle spins [60]. To calculate the probability of this chain 
existing in the LVS we find the fraction of points which have the mass ordering 
rriqj^ > m^o > mj^ > m^o. We find that this ordering occurs in all the points we 
generate, so that this chain should be useful for analysing LHC data. 
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4.2 The Dark Side 

We now discuss the possibility that /i < 0. It is well known that there is a correlation 
between the sign of /i and the sign of the SUSY contribution to the anomalous 
magnetic moment of the muon in constrained models like mSUGRA and the LVS 
(although not in the general MSSM). Hence while experiment currently favours fi > 
0, and it is true that a negative (5a^ leads to a large value, it could be possible 
to offset the bad fit to (5a^ with particularly good fits to other observables. Also of 
interest is the fact that it is possible to satisfy the B-term condition in eq. ( p.l4| ) 
for some limited regions of parameter space if /i < and ~ lO^^GeV. When 
this condition is satisfied the LSP is the stau, which is ruled out by anomalous 
isotope abundance constraints. However it is conceivable that when the flux terms 
are computed and taken into account that this might change, and it might be possible 
to satisfy all the equations in ( 2.14| ) in a phenomenologically viable way, without 
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introducing extra TeV scale matter as discussed above. We therefore press on in our 
exploration of the dark side of the Large Volume Scenario. 



(a) (b) 
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Figure 8: Posterior pdfs for /x < in the mo-tan /? plane for (a) flat priors and (b) REWSB 
priors with 68% and 95% c.l. contours shown. 

Fig. H shows the likelihood distributions marginalised to the mo-tan /3 plane 
for both sets of priors with 68 and 95% c.l. contours. The pdfs display a similar 
"boomerang" shape as Fig. ^ does in the /i > case, except the viable region of 
parameter space is considerably smaller. The posterior pdf for the fiat priors exhibits 
some slight bimodality, which is eliminated by the REWSB priors' pull to lower tan (3. 
The favoured region in both cases is around the centre of the boomerang, with the 
REWSB priors increasing the amount of likelihood at higher uiq. Similar to the 
yU > case, most of parameter space above and to the right of the favoured region is 
forbidden by having a stau LSP. 
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Table 7: Comparison of the best-fit points and statistical pulls for n < and ^ > with 
flat priors. niQ and tan /? are the values in parameter space, and all other values are x' 
values. The third column gives the error in our estimation of x'l ^ described Section ^ 
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Table |^ compares the best-fit points of botli signs of /i, for fiat priors. Tlie 
WMAP upper bound is saturated by tlie neutralino. Tlie liglitest Higgs mass falls 
just under the 114.4GeV 95% c.l. lower bound set by LEP2. {g — 2)^ is slightly 
negative, making a significant contribution to the overall x^- The 6-observables 
BR{b 57) and Aq- are fitted slightly better by /i < 0, while for AMb^, sinO^Q 
and M\Y there is essentially no difference. The branching ratio BR{B /U^/U^) 
comfortably evades experimental bounds in both cases. As in Section ^ we present 
the error in our estimation of the x^s of the observables, and the total error in y^. 

The overall of the fit is 20 and 20.29 for the fiat and REWSB priors respec- 
tively. In both cases nearly 65% of the comes from {g — 2)^. If we omit this 
observable, then the is 6.25 (7.69), which is roughly the same as for > with- 
out {g — 2)^, which has x^ = 8-10 (5.77). Therefore, if were to omit the anomalous 
magnetic moment from the likelihood both signs of n would be consistent with the 
data. Calculating the p-values associated with the total x^s (including {g — 2)^) we 
obtain Pf[gj^{fJ' < 0) = 0.010 and -P[^eWSb(/^ < 0) = 0.009. This is statistically 
significant down to the 2% level, and allows us to reject the possibility that /i < 0. 



5. Profile Likelihoods 

A regular objection to the Bayesian anal- 
ysis presented above comes in the form of 
criticism about the the subjectivity of the 
priors. While we have shown that in the 
LVS the posterior pdfs are essentially inde- 
pendent of the priors, we believe it fair to 
present a frequentist analysis of the model. 
Frequentists prefer to get rid of nuisance 
parameters by maximising them. This is 
known as concentration of parameters and 
the likelihood function of the reduced pa- 
rameter set is called the profile likelihood. 
However, it is important to state that pro- 
file likelihood plots are not pdfs, as the pro- 
file likelihood is not derived from a proba- 
bility distribution. We can easily derive the 
profile likelihood from the Markov chains we have generated as follows [8]: we bin 
the chains in the usual way, and then find the maximum likelihood in each bin and 
plot that. The 95% and 70% c.l. regions are then defined by 2AlnL = 5.99 (2.41) 
respectively where AlnL = In L^ax — InL [61]. 
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Figure 9: 2D profile likelihoods in the 
m-o-tan/? plane for ^ > 0, including 70% 
and 95% confidence limit contours. 
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In Fig. 1^ we plot the profile likelihood with 70% and 95% c.l. contours in the 
mo-tan /3 plane. The plot is reassuringly similar to those in Fig. ^. This illustrates 
an important point: that given enough data we expect the profile likelihood and 
the Bayesian likelihoods to look the same. Fig. ^ shows profile histograms for some 
relevant sparticles. The profile histograms have been rescaled to ease comparison 
with the Bayesian posteriors on the same plots. The profile results are in good 
agreement with the Bayesian histograms, further illustrating the point made above 
regarding convergence of profile and Bayesian likelihoods in the presence of ample 
data. The tails of the histograms are slightly noisier than those of the Bayesian 
posteriors, but this could be eliminated by running the MCMCs for longer. 

Using profile likelihoods we can establish the nature of the spike feature discussed 
in Sect. ^. We show in Fig. ^ the profile likelihoods for VtoKih'^ and for the mass 
splitting — m^o , which have been multiplied by constants for comparison with the 
Bayesian posterior distribution with flat priors. Fig. ^(a) shows that the region of 
very low relic density does not fit the data any better than other regions of parameter 
space. This establishes that the spike in the relic density is a volume effect, as 
promised earlier. In Fig. |lO](b) we see that the peak in the stau-neutralino mass 
splitting is also a volume effect, which was shown in Fig. ^to be associated with the 
spike. 




Figure 10: Histograms for (a) dark matter relic density and (b) — m^o for > 0. 
For tlie "flat" plot we show the posterior pdf per bin. For the "profile" plot we show the 
likelihood profile discussed in the text. To aid visual comparison the profile plots have 
been multiplied by constants. 

6. Conclusion 

We have used Monte Carlo Markov chain methods to make global fits to the Large 
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Volume Scenario in the minimal case with modular parameter A = 1/3, the first time 
this method has been applied to a model derived directly from string theory. As 
indirect constraints on the model we have used the WMAP 5-year dataset, the most 
recent measurement of rrit from the Tevatron and a suite of other electroweak and 
B-observables sensitive to flavour changing neutral currents, including for the first 
time in MCMC fits BR{B — > ru), AMb^ and the isospin asymmetry Ao_. We have 
shown that the model is constrained enough that the choice of flat or REWSB priors 
does not make a radical difference to the posterior pdf, unlike the CMSSM case. This 
illustrates the point that although the priors encode our uncertainty or a priori beliefs 
about a quantity, given enough data nature will speak for herself and the posterior 
pdfs will become essentially independent of the priors. Furthermore, the frequentist 
profile likelihoods we have presented have identified the same region of parameter 
space as the Bayesian hkelihood, indicating the robustness of our fits which again 
contrasts with the CMSSM case. We have constructed a new quantitative measure 
to determine which observables are constraining the form of the likelihood the most, 
based on examining the difference between posterior pdfs when we include and omit 
an observable from the construction of the likelihood. We find that {g — 2)^, the 
Higgs mass and the dark matter relic density have the greatest effect on the fits. 

We have also investigated both signs of /i, and by calculating P-values from the 
of the best-fit points in both cases we have rejected the possibility that < 0, 
while n > is consistent with our constraints. 

One of the main contributions of this paper is in forecasting what might be seen 
at the LHC. We find that all the points we generate allow the "golden channel" 
cascade Ql ^ X2 ^ ^ Xi- Use of this cascade should allow us to extract much 
information about masses and even spins of observed sparticles. Also, the 95% upper 
bounds on several important sparticle masses all fall below 1.2 TeV within reach of 
the LHC, implying early SUSY discovery at LHC. However the probability that 
rUr — m^o < lOGeV is 38% so that reconstructing the stau may prove difficult, and 
the mass of the lightest Higgs nih is constrained to be less than 120GeV indicating 
that the Higgs may not be found at LHC for a number of years. A "smoking gun" 
signature for discriminating the model against the CMSSM could be the ratio of 
gaugino masses 1.5 — 2 : 2 : 6 for the LVS, compared to 1 : 2 : 6 in the CMSSM, 
as first noticed in ref. [6]. In our model tan/3 is bounded from above at around 
tan /3 = 20, in contrast with the CMSSM where the region with the highest likelihood 
occurs where tan f3 > 50. While we must take into account the inevitable fuzziness of 
our bounds due the fluxes, a precise measurement of tan/3, although difficult, could 
therefore prove sufficient to reject the minimal LVS. 

Further work in this vein could include taking into account the modular parame- 
ter A as an extra variable in parameter space, thereby indirectly probing the geometry 
of the Calabi-Yau on which the model is compactified. In terms of connecting with 
the LHC, now that we have successfully sampled the likelihood distribution it is pos- 
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sible to construct correctly weighted samples of collider observables, as done using 
randomly generated points in [6]. To escape from the dilute flux approximation it is 
ultimately desirable to take into account analytically the fluxes on the Calabi-Yau, 
thereby allowing fully realistic behaviour at the string scale. This is unfortunately 
easier said than done, but the effects this would have on string phenomenology should 
not be underestimated. 
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